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We introduce a new 2A'^-paranietric family of maximally superintegrable systems in A^ dimensions, 
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space. These systems possess closed bounded orbits and integrals of motion which are polynomial 
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I. INTRODUCTION 

The aim of this paper is to introduce a new class of 
maximally superintegrable systems that are obtained as 
a symplectic reduction of the anisotropic harmonic oscil- 
lator. These systems depend on a set of N real and N 
integer parameters and possess integrals of motion poly- 
nomial in the momenta. The Hamiltonian defining this 
family is 
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We recall that in classical mechanics, superintegrable 
(also known as noncommutatively integrable [2J|) sys- 
tems are characterized by the fact that they possess 
more than TV integrals of motion functionally indepen- 
dent, globally defined in a 2A^-dimensional phase space. 
In particular, when the number of integrals is 27V— 1, the 
systems are said to be maximally superintegrable. The 
dynamics of these systems is particularly interesting: all 
bounded orbits are closed and periodic. This issue, for 
the spherically symmetric potentials, was first noticed 
by Bertrand [3|. The phase space topology is also very 
rich: it has the structure of a symplectic bifoliation, con- 
sisting of the usual Liouville- Arnold invariant fibration 
by Lagrangian tori and of a (coisotropic) polar foliation 



llj . Apart from the harmonic oscillator and the 



Kepler potential, many other potentials turn out to be 
superintegrable, like the Calogero-Moser potential, the 
Smorodinsky-Winternitz systems, the Euler top, etc. 

A considerable effort has recently been devoted to the 
search for superintegrable systems as well as to the study 
of the algebraic and analytic properties of these models. 
For a recent review of the topic, see [32 1. 

The notion of superintegrability possesses an interest- 
ing analog in quantum mechanics. Sommerfeld and Bohr 
were the first to notice that systems allowing separation 



of variables in more than one coordinate system may 
admit additional integrals of motion. Superintegrable 
systems show accidental degeneracy of the energy lev- 
els, which can be removed by taking into account the 
quantum numbers associated to the additional integrals 
of motion. One of the best examples of this phenomenon 
is provided by the Coulomb atom Il2| m , [l9|, which is 
superintegrable in N dimensions [7|, [28|. A systematic 
search for quantum mechanical potentials exhibiting the 
property of superintegrability was started in 14 1, |22l | 
and [35|. These models in many cases are also exactly 
solvable, i.e. they possess a spectrum generating alge- 
bra, which allows to compute the whole energy spectrum 
essentially by algebraic manipulations [SlJ. In classical 
mechanics, the multiseparability of the Hamilton-Jacobi 
equation implies that there should exist at least two dif- 
ferent sets of N quadratic integrals of motion in involu- 
tion. Reduction techniques both in classical and quan- 
tum mechanics are well-known (see, for instance, [J]). 
Essentially, the common idea of several of the existing 
approaches is to start from a free motion Hamiltonian 
defined in a suitable higher-dimensional space and to 
project it down into an appropriate subspace. In this way, 
one gets a reduced Hamiltonian that is no longer free: 
an integrable potential appears in the lower-dimensional 
space [26| . A different point of view, that we adopt here, 
is to start instead directly from a nontrivial (i.e. not free) 
dynamical system in a given phase space and to reduce 
it to a proper subspace, in such a way that the superin- 
tegrability of the considered system be inherited by the 
reduced one. 

In this work, we study the reduction of an anisotropic 
harmonic oscillator, defined in a 2A''-dimensional classi- 
cal configuration space. This system is maximally super- 
integrable. It is described by the Hamiltonian 
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We will prove that it can be suitably reduced to the sys- 
tem ([l}, and that this new system is still maximally su- 
perintegrable, with integrals of motion inherited from the 
system 1^. This goal is achieved under the assumption 
ni = ^2, ...,n2N-i = n2N- From a geometrical point of 
view, the approach we adopt reposes on the Marsden- 
Weinstein symplectic reduction scheme [l|, [21j, [23|. 
Given a symplectic manifold (M, $7), let Ki,...,Kk be 
k functions in involution: 



{K,,K^} = t,3 = l,...,k. 



(3) 



Assume also that dKi are independent at each point. 
Since the flows of the associated Hamiltonian vector fields 
Xki , ■■•, Xxk commute, they can be used to define a sym- 
plectic action of G = M'' on the manifold. Let J be the 
momentum map of this action, and /i be a regular value 
for J. Then we can conclude that P^ := J"^(/i)/G is still 
a symplectic manifold, of dimension dimM — 2k, called 
the reduced phase space. In our case, Ki,i ^ I, ...,N are 
components of the angular momentum, J = Ki x ... x Kj^ 
is the momentum map, G = 50(2) x SO{2) x ... x S0{2) 
{N times), and the reduced space is P^ ~ J^^ {^)/T^ , 
where T^ is the TV-dimensional torus and dimP^ = 2N. 
This procedure is a generalization of what in celestial 
mechanics, since the work of Jacobi, is called "elimina- 
tion of the nodes" (see [l|] , chapter IX for details) . The 
reduced Hamiltonian is reminiscent of the structure of 
the original Hamiltonian, defined in the 4A'^-dimensional 
phase space, but also possesses a Rosochatius-type term 
[29|, [la |. involving parameters ki corresponding to the 
variables that become ignorable, in addition to the har- 
monic part. Therefore, using the reduction procedure, 
we obtain the parametric family of Hamiltonian systems 
([T]), defined on a reduced phase P^. 

The transformations we consider, although very sim- 
ple, are non-trivial, since the reduced Hamiltonian is not 
shape-invariant. Nevertheless, since the reduced system 
turns out to be maximally superintegrable, bounded or- 
bits still remain closed in the reduced space. 

For N = 2 maximal superintegrability lj|, 35| and 
exact solvability 3l| of the system ^ was already es- 
tablished for ni = n2 = 1 and for tii ~ l,n2 ~ 2, fc2 = 0. 
The integrals of motion in these cases are second order 
in the momenta. 

Here we will show that in the general case (n^ and N 
arbitrary positive integers and ki arbitrary real numbers) 
N integrals can be chosen to be of order two, the other 
A^ — 1 functionally independent ones of order Ui + n^ or 
iT-i + UN — 1- Other systems possessing third and higher 
order integrals have been studied in the literature [Q|, 



m, [27 



This paper is directly related to the recent interesting 
work by Verrier and Evans 3J] , that performed a similar 
reducing transformation for the Kepler potential. They 
found a superintegrable system in three dimensions pos- 
sessing a quartic integral. They also conjectured that 



the system ([T]) in three space dimensions should be max- 
imally superintegrable, although the explicit expression 
of the integrals remained to be determined. In the fol- 
lowing, we will prove this conjecture, and also we will es- 
tablish that the system ^ is maximally superintegrable 
in full generality, i.e. for TV arbitrary, providing explicitly 
the corresponding set of integrals of the motion. 

After the present article was submitted we learned of 
a new article by Evans and Verrier [lO| in which the au- 
thors also establish the superintegrability of the system 
^ ioi N = 3. Their results are compatible with ours, 
though they express the integrals in terms of Chebyshev 
polynomials. Moreover, they also treat the quantum ana- 
logue of system ([1]) and establish the degeneracy of the 
energy levels related to the representation theory of the 
group SU{3). 

The paper is organized as follows. In Section II, the 
main properties of the anisotropic oscillator are briefly re- 
viewed. Then its reduction to the planar case is studied 
in detail. We will show how superintegrability is pre- 
served under a multipolar change of variables and sub- 
sequent reduction. In Section HI, the same problem is 
treated and solved in full generality. Some open prob- 
lems are discussed in the final Section. 



II. REDUCTION OF THE ANISOTROPIC 
OSCILLATOR 

The anisotropic oscillator in the two-dimensional case 
both in classical and quantum mechanics was discussed 
by Jauch and HiU ^, ^, ^. The system Q is also 
known to be superintegrable in any dimension, if the ra- 
tios of the frequencies are rational. Let us consider a 2N 
dimensional space and assume 
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Following [19|, we define the set of invariants in an 
auxiliary complex phase space, with coordinates Zi,Zi, 
i = 1, ..., 2N. Precisely, 

Zj = Pj ~ injujyj, Sj = pj + in-jUiyj. (5) 

It is easily checked that the expressions 



Cjk 






(6) 



provide integrals of motion. They can be also arranged in 
a real-valued form, as the combinations (1/2) (cy +Cij) 



and (l/2i) {a 



In particular, among these inte- 



grals we have the angular momenta 

Lik = Vipk - TJkPi 
(when Ui — rik) and the tensor 



T,,, 



PiPk + niUkUJ y^yk 



(7) 
(8) 



We will now study reductions of the anisotropic oscilla- 
tor ([2]) and establish the superintegrability of the corre- 
sponding dynamical systems. 



Hamiltonian and first integrals: the planar case 

We recall the definition of a momentum map. For fur- 
ther details, see for instance [l|]. Let (M, il) be a 2n di- 
mensional symplectic manifold. Suppose that a Lie group 
G acts on M and leaves Vt invariant. Let g be the Lie 
algebra of G, g* its dual space, and <,> the natural 
pairing between the two spaces. 

A momentum map for the G-action on (il/, f7) is a map 
J : M ^r Q* such that, for all X e g, 

d{{J,X))=ix^ 

In particular, if the manifold is exact, i.e. f2 — d0, and 
the G-action leaves 9 invariant as well, we have 

Jx = ixd 

We will also assume that the map is equivariant with 
respect to the coadjoint action Ad* of G on g*, i.e.: 

for aU gr G G, f e q* and X e q. 

Let us first consider a simple case, when the anisotropic 
oscillator is defined in a symplectic manifold M with 
dimM — 4. So, fi = Si=i dy"^ ^ dpi. In order to make 
the reduction possible, we select frequencies to be equal 
in pairs, so that we have only two independent frequen- 
cies. Hence the system ^ takes the special form 



Hi = 2(Pi+P2+l'3+P4) + — ^(^1+^2) + ^— (2/3+2/4)■ 
(9) 

In the auxiliary coordinates zi,Zi,... ,24,24, we have ex- 
plicitly 



21 = Pi -miujyi, Z2^P2- 


iniujy2, 




Z3 = P3 - i"2W2/3, 24 = P4 - 


[n2Ujy4. 


(10) 


luently, the Hamiltonian reads 
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Put in a matrix form, the set of invariants ^ can be 
represented by the matrix 
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Let us consider now the following change of coordi- 
nates: 



Xi cos 0:3, 

X2 COS 0:4, 



y2 = a;isma::3 
2/4 = X2 sin 0:4. 



The corresponding momenta read 
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(13) 



Pismxa, 

' P2 sin X4 . 
(14) 
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Xn 
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x\ 


=X4 + 02 



The group T2, which is the group 5*0(2) x 5*0(2) in 
the old coordinates, acts on R'' as follows 



(15) 



This group leaves fi invariant. The fundamental vector 
fields on r*R'' corresponding to this action are: 

Xi^d^,, X2^d,, (16) 

and, ii X = AiXi-}-A2A'2, the momentum map J satisfies: 

J{ai,a2) = ^(Al9x3 + MdxJ = A1P3 + A2P4 (17) 

Let us choose a regular point in t* (the dual of the Lie 
algebra of T2), for instance 

(18) 



= Vk 



P3 = Vki, Pi 
The inverse image under J is 

J^^{\/h,\/%) = {pi,P2, \fk[, \fk'2,Xx,X2.,X3,X4) 

(19) 
The stabilizer of this point in ^3 under the coadjoint 

action of T2 is the whole group, because its action is 

trivial on the p coordinates. 

The reduced phase space is therefore 

J~'(yfc^, yfc^)/T2 « {(pi,P2, xx,X2) e R4} (20) 
and the reduced Hamiltonian is 



pI , P2 . fcl , ^2 



7-7- r' I I ^z I '"-I- ( '"^ I i Z Z > z z Z /m\ 

-^2 - ^ + ^ + :73 + :73 + ^^ ^1 + ^^ 2^2 (21) 
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Let i^ be a first integral of the Hamiltonian H4{p,y), 
i.e. {H4,F} = 0. 

We show now how the original ring of integrals can 
be reduced in the low-dimensional phase space. First, 
we consider the restriction F of the function F to the 
manifold J~^{\/hi, v^)- 



Observe that F can be defined on the quotient mani- 
fold J~^{^/kl, y/k2)/T2, when it is constant on the equiv- 
alence classes, that is, -F is independent on X3,X4. In this 
case F can be factored out in the following way: 




where tt is the canonical projection and 



Fr0 7r = F. 



(22) 



Then, 



{H2,Fr}^0. (23) 

The integrals of the system (jlip are given in the matrix 
P^ (though only 7 of them can be functionally indepen- 
dent). Those that will survive the reduction (|^D|) are the 
ones that are left invariant by the S0{2) x 5*0(2) rota- 
tions (jl5|l . They must Poisson commute with 



Reduction of the first integrals 

The reduction is performed using the change of vari- 
ables (|13|) . (fT4|) and the convention (fT8|) . The integrals 
([M)l reduce to constants L12 = ^/ki, L34 = y/h^. The in- 
tegrals (|29|) reduce to nontrivial integrals for the Hamil- 
tonian in Eq. (I2ip . namely 
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(21Z2 - 2;22;i) = yip2 - y2Pi, 

(Z3Z4 - 24^3) = y3P4 - y4P3. (24) 



2n2a; 

The Poisson bracket can be written in terms of the z, 
variables as 



{f{zi,Zi),g{zi,Zi)} 
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N 2k 
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df dg df dg 



(25) 



, dzi dzi dzi dzi 

k=lj=2k-l \ J J J J, 

(in this Section we have iV = 2). 

Functions of Zk,Zk Poisson commuting with L12 and 
L34 must satisfy 

Z2dzif - zidz^I + Z2dzif - zidz^f = 0, 

Zidz^ f - z^dzj + Zidg^ f - zsdgj = 0. (26) 

A basis for the corresponding SO{2) x SO{2) invariants 
is given by 



zl, Ci=zf + zl, r]i = zizi + Z2Z2, 



^i=zl 

£,3=zl + zl, ^3 = z|-fZ4, 112 ^ Z3Z3 + Z4Z4, (27) 

Finally the integrals of motion must satisfy 

{i^4,/(Ci,6,r?i,e3,6,r?2)} = 0. (28) 

Solutions of equation ([28]) are for instance 

\2\ 



F, = \{\z,\' + \z2n, E2^li\z,\^ + \Z4\ 

Q, = {zi + zir{zi+zir\ 
Q, = izi + 4r{zi+zir\ 

h = {zl + z^){zl + zl), h = {zl + zl){zl + zl). 



(29) 



Only five of these integrals are functionally independent. 



The remaining two integrals in (|29p give nothing new and 
we have 

h=4:{Ef-kinluj^), l2 = 4:{Fi-k2nluj'^). (31) 

Three functionally independent real integrals of motion 
of the system with Hamiltonian pTj) can be chosen to be 



(32) 



{Fi,E2,Q = ^{Qi + Qi)}. 



They are of order 2, 2 and 2(ni + n2) in the momenta, 
respectively. Their existence is the proof of the maximal 
superintegrability of the considered system. 

The integral of motion Q simplifies to give a second 
order one in two cases (that were known previously [ij] , 
H). They are 



I) ni = ^2 = 1 
4S1S2 - Q 



2lu^ 

II) ni =: l,n2 = 2,/c2 
'8 FfF2~Q 
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The integrals ([33]) and ([34]) are responsible for the sep- 
aration of variables in polar and parabolic coordinates, 
respectively. The integrals {Fi,E2} are responsible for 
the separation in cartesian coordinates. 



III. THE GENERAL CASE 



where 



Within the same approach, it is easy to extend the 
previous picture to the general situation of a reduction 
from a 2N to a A^-dimensional configuration space: 



2N 
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Indeed, let us introduce the affine variables 

Zk = Pk - T^nkujyk, k = l,...,2N 
so that the Hamiltonian reads 
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The Poisson bracket is defined as in eq. (|25p . The invari- 
ants under the S0{2) x • • • x SO{2) group action gener- 
ated by Li2, • • • , L2N-1.2N are 



C2fc-1 



2 (2 

^2fe-l + ^2fc7 

-2 , -2 

^2fe-l + ^2fc 
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,iV, 



(36) 



apart from the quantities L12, • • • , L2n~i,2N and the "2- 
plane energies" which commute with the Hamiltonian 
H2N 



Imposing 



\zi\^ + \z2\'^,...,\z2N-i\^ + \z2n\^ 
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(38) 



where i = (Ci, ■ • ■ ,6Af-i), C = (Ci, • ■ ■ :6jv-i), we get 
the differential equation 
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Its general solution depends on 2N — 1 invariants, which 
can be chosen as: 
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(40) 



Q2fe-i and J= I^J + ^al- 

Using the transformation P^ we now reduce the orig- 
inal Hamiltonian to the following one: 



1 ^ 1 ^ A- UJ^ ^ 
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The corresponding reduced invariants are: 



1 



k 



1 



^' "2^' ^ 2;?+ 2''''^ "^^ 



^2/-i — :r (Q; + Qi) , 



Z = l, 






,N, 



(42) 



/ 2 , "-' 222 o- M^ 

/ = (P; + — ~ '^i^ 2;; - 2iniu>pixi) 



X (p' 



«JV 



"-AT 



tt-jyCj x^ -I- 2iniqujpNXNY 



There are 2Af — 1 functionally independent integrals and 
consequently the system is maximally superintegrable, 
proving the conjecture of |34l |. 

We have thus added a new maximally superintegrable 
system in A^ dimensions to previously known ones (see, 
e.g., i, i, 0, M. S). 



IV. OPEN PROBLEMS 

From the previous considerations, it emerges that it 
would be desirable to construct systematically transfor- 
mations mapping a superintegrable system into another 
system, that is also superintegrable, and defined in a re- 
duced phase space. It seems natural to associate such 
transformations to the rich symmetry structure possessed 
by superintegrable systems. For instance, changes of vari- 
ables of the type (fT3|) are clearly related to invariance 
properties under rotation. From this point of view, the 
role of higher order groups of transformations generated 
by the flow associated to integrals that are polynomi- 
als in the momenta remains to be fully investigated. A 
quantum mechanical version of this reduction procedure 
is also to be understood. For A^ = 3 the quantum system 
was treated in [1(|. The reduction was performed for the 
classical system. The reduced system was then quantized 
in cartesian coordinates. 
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